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Abstract 

We provide evidence of the relation between snpersymmetric gange theories and matrix 
models beyond the planar limit. We compnte gravitational conplings in gange theories 
perturbatively, by snmming genus one matrix model diagrams. These diagrams give the 
leading 1/A^ corrections in the large N limit of the matrix model and can be related to 
twist held correlators in a collective conformal held theory. In the case of softly broken 
SU{N) Af = 2 super Yang-Mills theories, we hnd that these exact solutions of the matrix 
models agree with results obtained by topological held theory methods. 


November, 2002 


1. Introduction 


Recently it has become clear that holomorphic or F-term information in A/” = 1 
supersymmetric gauge theories can be exactly computed using perturbation theory, when 
these terms are considered as a function of the glueball superheld S [Q (for earlier work see 
among others . Furthermore, for a large class of theories that allow a large N 

expansion a la’t Hooft 0 , the held theory Feynman diagrams, computed in a background 
super gauge held, can be seen to reduce directly to the diagrams of a zero-dimensional 
bosonic matrix model, where the matrix model potential is given by the gauge theory 
tree-level superpotential 0 ; for a recent alternative derivation of this fact using anomalies 
see 


10 


One such holomorphic quantity is the ehective superpotential lFeff(*S') that is given 
by summing just the planar diagrams, even for a hnite rank gauge theory. Non-planar 
diagrams will in general contribute to gravitational corrections [|^. In particular, diagrams 
with genus one topology, that give the leading 1/A^^ correction J^i to the matrix model 
free energy, contribute to an ehective curvature term of the form 
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IGtt^ 


d'^x J^i{S)TtR+ AR+, 


( 1 . 1 ) 


with i?_(_ the self-dual part of the Riemann curvature tensor. This induced gravitational 
correction measures the back-reaction of the held theory when it is placed in a curved 
background. The actual derivation along the lines of [P] of the matrix model Feynman 
rules, directly from the gauge theory Lagrangian in a curved superspace, will be presented 
elsewhere 0 - 

For exactly solvable matrix models the summation of the diagrams, of any hxed topol¬ 
ogy, in closed form can be done in principle, although the techniques become progressively 
cumbersome for high genus. One can thus try to compare these exact answers to known 
properties of four-dimensional supersymmetric gauge theories. 

In this note we will compare the results for a single matrix model to the gravitational 
corrections that have been computed for topological held theories that are twisted versions 
of A/” = 2 supersymmetric Yang-Mills theories |0JT3,0,|T^ . These topological held theories 


are used to compute the Donaldson and Seiberg-Witten invariants of four-manifolds. The 
gauge theoretic results have been derived making use of the Seiberg-Witten solution 


and holomorphy and duality arguments. In this paper we will demonstrate how these 
terms can also be computed using loop equations of matrix models. 
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The genus one correction is simpler for many reasons. Being a first order correction 
to a semi-classical evaluation it is given by a fluctuation determinant, also in the matrix 
model. If the matrix model is exactly solvable, this gets often reflected in an emergent 
geometry, that in terms of topological string theory arises from a geometrical transition 
from an open string to a closed string description. General arguments tell us that such a 
dual geometry takes the form of a non-compact Calabi-Yau three-fold. Topological closed 
strings propagating on such a CY three-fold give rise to a genus one partition function T\ 


that can be expressed as a generalized Ray-Singer analytic torsion |T7 


Tx = pq\ 

p,q=Q 


1)P+'? logdet A 




( 1 . 2 ) 


where Ap^q = {9, 9 } is the Laplacian acting on (p, q) forms. In the simple class of matrix 
models that we consider in this paper the effective geometry is essentially given by an affine 
algebraic curve, and therefore we expect for general reasons an expression of the form 

J’l = -^ logdet Ao 

with Aq the (scalar) Laplacian on the algebraic curve acting on the collective bosonic held. 
We will verify this is indeed the case in some cases by explicit computation. This relation 
between matrix models and two-dimensional conformal collective held theory is a much 


more general feature, see e.g. [18 


The plan of this paper is the following. In section 2 we state the precise relation 
between the matrix model and the gauge theory quantities relevant for gravitational cor¬ 
rections. In section 3 we consider gravitational couplings in A/” = 2 super Yang-Mills 
theories obtained by a topological twist, and make a comparison to the matrix model re¬ 
sult for SU{2). Then in section 4 we analyze the general matrix model answer in terms 
of a collective conformal held theory. We hnd that the genus one contribution can be de¬ 
scribed in terms of twist held correlation functions with extra dressing to match the loop 
equations. This will allow us to make a precise identihcation for general SU{N) gauge 
group. 

After this note was hnished appeared that discusses similar issues. 
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2. Superpotentials and gravitational couplings 

Let us first briefly review the main results of [P,|l[] for the prototypical case of a U{N) 
gauge theory coupled to single chiral matter held $ in the adjoint representation. We start 
with a tree-level superpotential 

j d^xcfe TrW{^), ( 2 . 1 ) 

where the polynomial W has n distinct critical points. If we consider a classical vacuum 
where one distributes W of the eigenvalues of $ in the fth critical point of W, we have a 
classical breaking pattern 


U{N)^U{Ni)x---xU{Nr,). 


The strong coupling dynamics of the corresponding quantum vacuum is captured by the 
effective superpotential Weff(5'i) as a function of the glueball superhelds 
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327r2 


Trsc/(Ari)>V^. 


According to the prescription of [|T| this effective superpotential is given by 

dTo 


WMS) = 


Ni 


dS, 


+ 27rro Si 


( 2 . 2 ) 


where tq is the bare coupling and J^o{Si) is the free energy of the corresponding matrix 
model, obtained in a semi-classical expansion around the classical vacuum. 

This matrix model takes the form of an integral over a N x N matrix $ (here we 
carefully distinguish between N, the rank of the matrix model, and iV, the rank of the 
gauge theory) 


-exp 

TrW($) 



= exp 


volf7(A^) J 

L 9s J 


g>0 


with the identihcation Si = gsNi in the ’t Hooft limit gg 0, Ni ^ oo. More precisely, 
we have 


rois) = ^sHog{s/A^„)+rr‘{s)- 


(2.4) 


The hrst term gives rise to the Veneziano-Yankielowicz effective action of the pure Yang- 
Mills theory [^, 

Wes{S) = NSlog{S/A^). (2.5) 
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Fig. 1: The simplest genus one, non-planar diagram in a cubic theory — the leading 
perturbative contribution to T\. The •’« indicate insertions of the background 
gauge held Wa. 


In the matrix model this contribution to (p.4|) is reproduced as the large N volume of the 
unitary group ||]. 


The second term in (|2.4|) is given by a sum over planar diagrams that appear in the 
perturbative expansion of the matrix model. (See m for a careful description of this 
expansion around a vacuum with a spontaneous broken gauge symmetry.) A diagram 
with £ index loops comes with a factor of S^. The actual physical values of WeS and the 
condensates Si in the quantum vacua are given by extremizing (|2.2|) with respect to the 
glueball helds Si. 

As we mentioned in the introduction this relation is not restricted to planar diagrams. 
There is an elegant interpretation of the higher genus diagrams that give the corrections 
jFg’s in terms of the coupling to a supergravity background [|l|. In particular, the induced 
gravitational effective action obtained by putting the held theory on a curved space-time 
contains the F-term 

^ j J^i{S)TtR+ a R+. (2.6) 


where Ry is the self-dual part of the Riemann tensor. (There is of course a similar anti- 
holomorphic term Ri multiplying TrR_ A R-.) 

If we consider the partition function on a Euclidean four-manifold then this 
gravitational coupling induces a term 


expri(S)(^x 



(2.7) 


with X the Euler number and a the Hirzebruch signature of M. 

Evaluating the term .T’i(5') in perturbation theory, one hnds that it is given exactly 
by the sum of diagrams with topology genus one, i.e. the diagrams that give the leading 
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l/N"^ corrections in the large N limit of the matrix model. More precisely, is given 

as 


+.^r‘(S)- (2.8) 

i 


This expression is the gravitational analogue of (O). Assuming conhnement, so that the 
only the held accounted for is S, the hrst term has an interpretation as an integrated form 
of the gravitational contribution to the 17(1) R-anomaly, 


8 ^ 

~ 167r2 


-Tr 

2 


U{N., 


aF AF - —TrRA R 
12 


(2.9) 


If one assumes that the low energy dynamics of the gauge system is described by an 
effective action in which the glueball superhelds Si can be treated as elementary helds, the 
anomalous behaviour under the R-symmetry 


S 




is reproduced by the combination of the Veneziano-Yankielowicz contribution A^5'log5' 
to WefF(5') — recall that the top component of S' is Tr F+AF+ — together with the 
— ^ log S multiplying the gravitational correction. Including the complex conjugated term 
that muliplies logS, we see that we pick up precisely the anomaly (|2.9|). 

The perturbative contribution to (^) is given by summing all genus one dia¬ 


grams. For example, in a cubic theory, with superpotential W ($) = the leading 

diagram is given by hg. 1 and this gives 


2 


Of course in the physical vacua all these expressions for Fi have to be evaluated 
for those values of the Si that minimize the effective superpotential given by the planar 
contribution. 0 

^ The genus zero diagrams can also contribute, but their contribution can be shown to cancel 
at the critical point. 
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3. Matrix models and J\f = 2 theories 


Matrix model methods can be used in particular to find the celebrated solution to the 
pure M = 2 super Yang-Mills theory of Seiberg and Witten . To describe the SU(N) 
gauge theory one breaks the supersymmetry down to Af = 1 by introducing a degree iV -|-1 
tree-level superpotential Tr W ($) of the adjoint chiral multiplet $ and picks the breaking 
pattern 

U{N) 


by distributing the N eigenvalues of $ equally among the N critical points of kb; that is, 
one chooses all Ni = 1. One further decouples the diagonal U{1) by putting the overall 
bare coupling tq = 0. The effective superpotential then simplifies to 

i 

In this case the planar diagrams can be exactly summed and the solution can be 
written in terms of period integrals on the associated hyperelliptic Riemann surface |^] 

N 

= P{xY + f{x), P{x) = W'{x) = UiX^~\ 

i=Q 

In this case the definition of the variables Si,..., Sn is subtle, since they are defined in 
terms of the traceless piece of a U{1) gauge field. Classically they vanish but as operators 
they make sense quantum mechanically [^. The dependence on the Si is implicit in 
terms of the quantum deformation /(x), a polynomial of degree N — 1. After solving the 
constraint dWeff{S) = 0, this curve takes the familiar SW form 


y2 = P(x)2 - A2^. 


(3.1) 


By the introduction of the bare superpotential we have effectively localized to a particular 
point of the Coulomb branch of the Af = 2 theory. 

To obtain the original Af = 2 model one can now scale the tree-level superpotential 
as kb —e kb and take the limit e —0. There are two obvious quantities that by a scaling 
argument do not depend on e and can therefore be straightforwardly extracted from the 
Af = 1 solutionil. First, there is the coupling matrix 

d^Po 

dS,dSj 

^ For the cubic potential kb = this can be explicitly checked by scaling m —> em, 

g ^ eg, S ^ eS. Note that in the matrix model S = gsN and therefore gs —> egs- 
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of the U{1)^ low-energy effective Abelian theory. Geometrically this is given by the period 
matrix of the cnrve (|3.1| ). The second e-invariant qnantity is the genns one free energy jFi 
that gives the gravitational correction (p.6| ). 

3.1. Gravitational coupling from topological field theory 

On flat spacetime the Af = 2 SU (2) gauge theory is described by the Seiberg-Witten 
solution. Putting the theory on a curved manifold additional gravitational terms appear in 
the low energy effective action. This gravitational correction has been directly computed 
in the M = 2 theory — more precisely, in a topological twisted version of the theory 
that computes Donaldson invariants. In the twisted version one modihes the action of the 
Lorentz group 

50(4) ^ SU{2)+ X 50(2)_. 

One replaces SU{2)j^ with the diagonal subgroup of 5C/(2)_|_ x SU{2)ii, where the last 
factor is the M = 2 internal R-symmetry group . 

In the twisted topological theory considered on a curved four-manifold M these inter¬ 
actions are restricted to the topological R A R* and R A R terms (with R± = ^{R± R*)) 
proportional to the Euler number x(M) and the Hirzebruch signature a{M) respectively. 
The gravitational couplings contribute to the partition function with the factor 


exp lb(u)x + c(u)a ], 


(3.2) 


where b{u) and c{u) are functions of the parameter u on the gauge theory moduli space. 

The precise form of the functions 6('u), c{u) can be inferred from analyzing the modular 
transformation properties of the quantum theory on the curved manifold. Cancellation of 
the modular anomaly and additional input from the singularity structure of Seiberg-Witten 
moduli space determines the measure contribution 

To co nn ect these computations in topological held theory to the physical theory we 
recall that for manifolds with metrics of SU{2) holonomy (hyper-Kahler manifolds) the 
topological twist is invisible since there is no holonomy in 517(2)_|_. We can therefore 
directly compare to the physical gauge theory. In that case the metric is pure self-dual, 
and we have 


a = -3X. 
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For example one could take M = K3 for which x = 24 and a = —16. So the overall 
contribution to the path-integral is 


exp 


(6(n) - -c{u))x 


If we compare this to (12.71) , where we use that for a self-dual geometry ~ fc’’ = X; we 
have the following identihcation between the matrix model and gauge theory quantities 

= h{u) - ^c(n). 


We will now check this relation in a number of cases. We will for convenience put x — 1- 


3.2. The M = 2 SU{2) theory 

In this case the Seiberg-Witten geometry can be described by deforming the M = 2 
theory with a tree level superpotential, 

W'($) = e($2(3.3) 


(For more details about the perturbative derivation of this particular case see m-) As 
described above, extremization of the effective glueball superpotential gives the Seiberg- 
Witten curve for SU (2) 

= {x^ — u)‘^ — 1. (3.4) 


Here the scale A is set to one for convenience and the factor e is absorbed. As we mentioned 
the physical quantities Ti and the coupling matrix rij are independent of the deformation 
parameter e. The curve has four branch points at 


Xi — Aa/w ± 1- 


(3.5) 


It is described by the two-cut solution of the matrix model with the potential IF($) given 
by (H). 


The genus one free energy for two-cut solutions in matrix models have been explicitly 
computed. Here we use the relevant solution of Akemann |^, which is an elaboration of 
the methods of [E4|, 


•^1 - ^log|^(^)l 

i=l 

- Xa + I log|(li - X2,)(X2 - Xi)\. 


(3.6) 















This solution was derived by an iterative genus expansion of the loop equation; we discuss 
this further in the next section. Here A is the discriminant of the elliptic curve (|3.4|) 


A = = 64(n^ - 1), (3.7) 

i<j 

and K{k) is the complete elliptic integral, where the nome k is expressed in the modulus 
r of the SW curve. The solution also depends on the hrst moments of the potential that 
are generally dehned as 



For the simple potential (|3.3| ) the contour can be deformed to inhnity, and one gets Mi = e. 

For comparison with the gauge theory result it is useful to express in terms of the 
parameter of SU (2) moduli space u. The elliptic parametrization of the SW curve (|^ ) 
can be written in terms of Jacobi 9 functions as 

I oA oS 

^ ^ 2+^3 -I = _ ^4 _ 

2 ( 02 ^ 3 )'’ 4 ( 02 ^ 3 ) 4 ’ 

where the dehnition of the 6 functions is as usual 

62 = ^ 

n6Z 

n6Z 

9^ = 

66Z 

with q = A useful identity they satisfy is 9\^9\ = 9^. The complete elliptic integral 

K{k) can also be expressed in 9 functions as 


K{k) = \9l 


With these elliptic parametrization the matrix model answer for the two-cut solution to 
T\ can be written as 


Tx 


1 1 4 

-logeT - log—— 

6 4 7r(e'2c'3) 


1 , 166 '! 

12 (^3)"' 


(3.9) 


The factor log e can be absorbed in the measure. 
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3.3. Comparison to the gauge theory 


We now have to compare this result to the topological held theory answer that reads 

iga 



= P(u^ - I)'/®, 


where a and (3 are constant coefficients. This contribution to the partition function should 
match with the matrix theory computation for the corresponding genus one contribution. 
To check this, it helps to rewrite the gauge theory contribution as 


^gauge 




(3.11) 


with 


da 

du' 


A = 64(^2 - 1). 


Here 





i f du 
drr \ da 


2 


is rewritten in terms of the “electric” period of the Seiberg-Witten curve a. Substituting 
the modular parametrization of the curve in terms of 6 functions we hnd 


A= 


(Ma)®’ 

Comparing with the matrix theory contribution 


du 2 ^ 


we hnd perfect agreement. 


3 . 4 . SU{N) generalization 

The gauge theory computation for the partition function can be generalized for the 
SU{N) theory. The generalization is based on a similar analysis of anomalies as for the 
SU{2) case. 

At a generic point on the Coulomb branch, where the gauge symmetry is broken to 
^(l)Ar-i, SU{N) theory can be described by the hyperelliptic curve 

2N 

y'^ = P{xf - 1 = - Xi), 

i=l 

N N 

P{x) = ''^^UiX^~'^, = 0. 

i=0 i=l 


10 



Here the are the symmetric polynomials of the roots of P{x), and Xi are the branch 

points of the curve. The hyperelliptic curve is a Riemann surface of genus g = N — 1. 

For a genus g Riemann surface one takes a basis of 2g homology cycles Bi) with 
canonical intersection product. The periods of the curve are then related to a set of dual 
holomorphic one-forms uji — x^~^dx/y as (we choose the homology basis slightly different 
then in [p^p!5| in order to make contact with the matrix model basis) 


A 





B. 


i-j 



daD,i 
duj+i' 


The period (or coupling) matrix is given as 


daD,i 


(BA-^),,. 


The partition function for SU{N) is a direct generalization of the corresponding SU (2) 
|. We will write it as (discarding overall constants) 


contribution [14 


^gauge 


— ^-x/2^cr/8^ 


(3.12) 


H = detHij, /S. = '^^{xi — XjY. 

i<j 

Putting a self-dual metric and y = 1 we get 

^gauge-A-l/2A-l/12. (3.13) 

To compare this result to the genus one free energy of the matrix model we hrst have to 
explain how matrix model results can be computed using conformal held theory. 


4. Multicut solutions and conformal field theory 

4^.1. Loop equations and Virasoro constraints 

For the one-loop free-energy for the SU{N) theory we have to solve the corresponding 
matrix model with the tree-level superpotential kF($) with kF'($) = eP($) and with the 
maximum number of cuts. The most efficient way to derive multicut solutions for matrix 
models is by using loop equations and conformal held theory techniques. The method of 
using loop equations to obtain the 1/N corrections in matrix models was developed in 
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Here we will follow closely that gives a good general exposition of the relation of these 
methods to conformal held theory. 

We start from the partition function of the associated matrix model 


Z = 


volU{N) 


exp 


TrW($) 

9s 


for the N X N matrix $ with a general potential W($). The reparametrization invariance 
of this integral leads directly to the so-called loop equationl. The simplest way to derive 
the loop equation is taking a shift $ 
equation for the loop correlator 


where h is a small number. This gives the 


where 


u(x)‘^ -^Tr 

9s 


u{z) = Tr 


X — $ 


N 


= 0 , 




X — $ X — Aj 


is the loop operator and A/ the eigenvalues of the matrix $. For a general potential with 
coupling constants 

OO 


n=l 


the loop equation can be rewritten as 

dx' 1 


j {r(x')) = 0. 


27rz X — x' 

where the contour includes all eigenvalues Aj but excludes the point x. Here we introduce 
the stress-tensor T(x) of the collective held (p{x) 


T(x) = ^{dip{x)f, 
ip{x) = W{x) — 2gs Tr log 


1 


X — $ 


The loop equation can be reformulated as the Virasoro constraints |25,26 


LnZ = 0, n > —1, 

Tl r\ OO 

^ o o , o 


dtk dtn-k 

k=0 " k=0 


dt 


n-\-k 


These loop equations have been recently given a gange theoretic interpretation in |M]. 
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where the operators satisfy (half of) the Virasoro algebra 


\Lin-i -bm] ^)-bTn+n- 

The loop equation can be solved iteratively order by order making a 1/A^ expansion. The 
planar limit is usually not so hard to solve, but for the next order solution one needs special 
techniques. 

4.-2. Planar solution and effective geometry 

To solve the loop equation in the N ^ (yo limit, it is simplest to rewrite it in terms of 
the polynomial 


If W(x) is of degree n + 1 then the polynomial f{x) is of degree n — 1. Denoting the 
classical average of the loop operator as 

1 

Ucix) = —{u{x)), 

N 

in the large N limit the loop equation becomes quadratic 


U)c{x)‘^ - —W {x)uJc{x) + 


Ag^N- 


-f{x) = 0. 


In terms of the bosonic collective held p{x) we can say that this obtains a large vac¬ 
uum expectation value in the planar limit, which is the solution for the classical Virasoro 
constraint. The current d(p{x) takes its classical value 

dipc{x) = W'{x) - 2gsNuc{x) 

= W'{x)-^^ + o(^) . 

X \x^ J 

For our solution we have with y = d(~pc{x) 


y‘^ = W\xf + f{x). 


We can write this as the hyperelliptic curve 


- Xi). 
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Since the classical field (f{x) changes sign around the branch points, the expectation values 
of the bosonic held are given by two branches of y{x). Then (p{x) can be thought of as a 


single bosonic held dehned on the branched covering given by the hyperelliptic curve (|4.1|) . 
The classical value of the collective held is given as 


dipc{x) = W{x - Xi)^/'^. (4.2) 

i 

4-3. Subleading corrections and twist fields 

The subleading term in the free energy is given by the Gaussian huctuations 
around the classical solution. As we mentioned in the introduction, by general arguments 
this term is equal to — |logdetAo, the Laplace operator on the Riemann surface, and 
additional (dressing) terms arising from the huctuation of the branch points. 

Instead of thinking of (p{x) as a held living on the hyperelliptic Riemann surface ([4.1|) , 
we can also think of it on the complex x-plane in the presence of twist operators a{xi) 
associated with the branch points Xi. In the neighbourhood of such a twist held the current 
dip{x) is no longer single-valued but has a branch cut in its operator product 

9(p(x) ■ a{xi) ~ (x - Xi)“^/^r(xi). 


Here cr(x) and r(x) are conformal helds of dimension 1/16 and 9/16 respectively. A naive 
expression for the genus one contribution to the matrix model would now be given by 

2N 

i=l 

The chiralS twist held correlation function is well-known iEzi@a 

ZtwiBt = .4^1 (4.3) 

i<j 


Here 


A = det(Aij) 


is the determinant of the period matrix, related to the integral of the one-forms uji over 
the A-cycles 


Aij — 


X 


3 ^dx 


(4.4) 


In this expression the chiral projection is done by putting the loop momenta of the held ip{x) 
to zero. 
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For example, in the case of a two-point fnnction we get the familiar resnlt 


{a{xi)a{x2)) = {xi - X 2 ) 


expressing the fact that the conformal dimension of a Z 2 twist held is 1/16. 

Formnla (|4.3|) also can be expressed as the chiral determinant of the Laplace operator 
Ao of the twisted boson on the hyperelliptic cnrve 


■^twist (det Aq^ 


-1/2 


4-4■ Star operators 


However, (|4.3|) is not the fnll answer, since as it stands this expression does not solve 
the Virasoro constraints. An elegant solntion to this has been given by Kostov in terms of 


star operators |]I8 


We can associate a Hilbert space with the local complex variable near each branch 
point and solve the Virasoro constraint in the vicinity of the branch point. We have to 
look for an operator which creates a conformally invariant state near the branch point. 
The twist operator itself does not satisfy all the Virasoro constraints, in particnlar it does 
not satisfy L_i. Therefore we will look for a new operator which satishes all constraints. 
Snch operators are called star operators [0, and they are constrncted from the modes of 
the twisted bosonic held near the branch pointi. 

The twisted bosonic cnrrent near the branch point Xi is now decomposed into a clas¬ 
sical and qnantnm part 


d(p{x) = dipc{x) -|- ar{x — Xi) 


—r—1 


rez+i 


The expansion of the classical cnrrent ( [1.2|) is 


\r—1 


dipc{x) = ^ HriXk) ■ (X -Xk) 

r>i 
' — 2 

This dehnes the coefficients The Fock vacuum for such a twist held is dehned as 

|0i) = (T(Xi)|0), 


^ We would like to thank I. Kostov for sharing with ns some nnpublished work on the con¬ 
struction of the star operators. 
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and it satisfies 


CKrIOi) =0, r > 0. 

Since it depends on the position of the branch points it is not translationally invariant. To 
make it invariant, one introduces the star operator 


S(xi) 






and assume it is dehned perturbatively by a mode expansion 

s(Xi) ^ ^ y ^ ^ri...rn ■ 

n>0 ' ri...r„ 

The coefficients in the mode expansion are determined by imposing the conditions of 
conformaf invariance 

=0, n>-l. 

Up to 1/N^ correction one hnds simply an extra multiplicative factor 

S(xi) = [n3/2ixi)] a[xi). 

The full genus one contribution to the free energy, obtained by solving the loop equa¬ 
tion including the order 1 /corrections, is therefore given by the correlation function of 
star operators, not the twist operators, 

2N 

= log( ) 

(4.5) 


= -^'^^Og^3/2{Xi) + log Zi 


i=l 

2N 


24 


twist 5 


i=l 


where -Z^twist is the correlation function of the 2N twist helds (|4.3|) . From the expansion of 
the classical current d(fc{z) we get 

1 1 

- -—logfliXi - Xj) 
i=l i<.j 

So the hnal result for ^matrix is theni 

^matrix = (4.6) 

This result is in complete agreement with the gauge theory partition function -^gauge (|3.13|) . 
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Note that this form is given incorrectly in |U]. 
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